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Abstract 

We continue our recently proposed holographic description of single-particle correla- 
tion functions for four-dimensional chiral fermions with Lifshitz scaling at zero chemical 
potential, paying particular attention to the dynamical exponent z = 2. We present 
new results for the spectral densities and dispersion relations at non-zero momenta and 
temperature. In contrast to the relativistic case with z = 1, we find the existence of 
a quantum phase transition from a non-Fermi liquid into a Fermi liquid in which two 
Fermi surfaces spontaneously form, even at zero chemical potential. Our findings show 
that the boundary system behaves much like an undoped Weyl semimetal. 
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1 Introduction 

Over the past years, the AdS/CFT correspondence has become a more and more popular 
and widespread tool which offers the opportunity to apply ideas from string theory to 
realistic materials studied in condensed- matter physics, see e.g. [1, 2, 3, 4, 5] and refer- 
ences therein. More specifically, it is potentially very useful when considering strongly 
coupled or critical condensed matter, which is generically not possible to describe within 
perturbation theory. However, in the context of the AdS/CFT correspondence these 
systems can be understood and investigated more conveniently using a dual theory in a 
higher- dimensional curved spacetime, i.e., in the framework of general relativity. Most 
commonly, the duality exists between a weakly coupled, classical gravity theory in a 
curved anti-de-Sitter (AdS) spacetime on the one hand, and a strongly coupled, confor- 
mal field theory (CFT) living on the flat boundary of the AdS spacetime on the other 
hand. Such a conformal field theory describes for example a quantum critical point of 
the condensed-matter system under consideration, and the observables that are most 
readily available from the AdS/CFT correspondence are usually correlation functions of 
the composite operators that classify the conformal field theory. For our purposes, most 
relevant are correlators of fermionic operators, since they may show Fermi or non-Fermi 
liquid- like behavior most easily [6, 7, 8, 9, 10]. 

However, in condensed-matter physics it is more natural to think in terms of fermionic 
single-particle operators, i.e., creation and annihilation operators in a Fock space, that 
satisfy the following (equal-time) anti-commutation relations 



where a labels the spin of the electron. Furthermore, in experiments on condensed- 
matter systems, the quantity that is measured is essentially always related to single- 
particle or two-particle correlation functions and not to the above-mentioned composite 
operators. For example, the retarded single-particle correlation or Green's function 
Gr is measurable in electronic systems using angle-resolved photoemission spectroscopy 
(ARPES). For this reason, we are interested in finding Gr of a strongly interacting 
condensed-matter system using holographic methods. Inspired by previous work on 
this topic, such as [11] and [12, 13], we have shown recently how the usual holographic 
prescription can be modified in such a manner that it allows for the construction of the 
retarded single-particle Green's function [14]. In a sense, our method can be seen as a 
bulk derivation of the semi-holographic description advocated in [12]. 



[lp a (x,t),t/jl l (x,t)] + = 5(X- x)5 a:a ' , 



(1.1) 
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A crucial consequence of the anti-commutation relations in equation (1.1) is the ex- 
istence of a sum rule for the corresponding single-particle two-point correlation function, 



This sum rule is essential in determining whether the quantity under consideration is 
the correlator of a single-particle or a composite operator, and therefore plays a central 
role in our construction. 

The quantum critical points described by the AdS / CFT correspondence are charac- 
terized by their dynamical scaling exponent z. In the usual correspondence, an anti-de- 
Sitter background leads to a quantum critical point with relativistic scaling, i.e., z = 1. 
But from a condensed-matter perspective we are also interested in quantum critical 
points that have a dynamical scaling exponent different from one. To achieve this, the 
usual anti-de-Sitter background is generalized into a so-called Lifshitz background, which 
leads to a non-relativistic, i.e., z ^ 1 scaling on the boundary [15]. Non-zero temperature 
effects can then be studied by placing a black brane in the Lifshitz spacetime [16], see 
also [17] for a more recent discussion. The pure Lifshitz geometry without a black brane 
develops singularities due to diverging tidal forces [18]; as a result the zero-temperature 
limit might become ill-defined 1 . While this could be of possible concern for our analysis, 
our findings show that the fermionic Green's function is well defined at zero temperature 
and the limit T — > is smooth, so at least the spectral-weight function for fermions does 
not suffer from any singularities. 

In [14] we have described how to construct the retarded Green's function of a strongly 
interacting, but particle-hole symmetric system of chiral fermions with an arbitrary dy- 
namical exponent z using a Lifshitz black-brane background. The aim of the present 
paper is to analyze the physics that follows from our prescription. In particular, we cal- 
culate the retarded single-particle Green's function in various cases, at zero and non-zero 
temperatures, and present the corresponding spectral-weight function and dispersion re- 
lations. Although we give some results for z = 1, we mostly focus on the case z = 2 
and four spacetime dimensions. With this number of dimensions, the chiral fermions are 
Weyl fermions and the boundary system turns out to be an interacting Weyl semimetal 2 
with Lifshitz scaling. By computing the momentum distribution of the particles and 
holes, we find for z / 1 that our holographic Weyl semimetal contains a quantum phase 
transition between a non-Fermi-liquid phase and a Fermi-liquid phase with two Fermi 
surfaces, one for the particles and one for the holes, even at zero doping. The scale set 
by the Fermi surface is determined by an appropriate combination of the dimension- 
less parameters A and g that enter in our prescription for the Green's function on the 

1 See however [19, 20] for possible resolutions within string theory. 

2 A semimetal is a gapless semiconductor. In addition, a Weyl semimetal [38] is a semimetal with 
touching valence and conduction bands based on chiral two-component fermions that in the non- 
interacting limit satisfy the Weyl equation ±<r • kip = Etp. Here, the ± denotes the chirality of the 



given by 




+oo 



(1.2) 



fermion. 
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(a) A > (b) A < 



Figure 1: An idealized picture of the ground state that we find for A > and A < 0. Grey areas 
refer to filled energy levels. In this picture, we only show the effect of the interactions on the 
occupation numbers and not on the dispersion relations. 

holographic boundary [14], and the quantum critical point is at A = 0. One can think 
of 1/A as a spin-orbit coupling constant and of g as an effective interaction parameter 
coupling the elementary fermion to the conformal field theory. The phase transition is 
schematically illustrated in figure 1, where we show the non-interacting band structure 
and how it is populated in the ground state of each phase. In particular, for A > the 
conduction band is empty and the valance band is completely filled, whereas for A < 
the conduction band contains a Fermi sea of particles and the valence band a Fermi sea 
of holes. Note that we present here an idealized picture that neglects the renormalization 
of the Fermi surfaces involved, which is discussed at length in this paper. 

The layout of the paper is as follows. In section 2, we briefly summarize some basic 
formulas from [14] that are needed here. In section 3, we present analytic and numerical 
results for the single-particle spectral function for the cases z = 1 and z = 2. In section 
4 we consider the single-particle momentum distribution, and describe properties of the 
quantum phase transition and the Fermi liquid involved, that does allow for well-defined 
but strongly renormalized quasi-particles and quasi- holes near the two Fermi surfaces. 
The main part of the paper ends with conclusions in section 5. Finally, there are a 
number of appendices which contain our conventions and give more details about the 
calculations. 

2 Single-particle Green's function for Lifshitz fermions 

In this section, we determine solutions to the Dirac equation in a Lifshitz spacetime, from 
which we construct the retarded fermionic single-particle Green's function. To keep the 
paper self-contained, we start by defining our notation and conventions. Readers who 
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are only interested in the result of the calculation may skip directly to equation (2.11). 
Further details on the derivation can be found in the appendices, in particular, for 
remarks on how to switch from natural to SI units we would like to refer to appendix 
A.l. Some of the results in this section have already appeared in more detail in our 
previous work [14] , other related work on fermions in Lifshitz backgrounds can be found 
in [21, 22, 23]. 

We begin by specifying the gravitational background metric of a Lifshitz black brane 
[16], i.e., 

ds 2 = _y2 {r) r 2 Zdt 2 + dr^ + r2df2 y2 =l _( r A\ d+z ~ l (2 .i) 

r A V z {r) \ r / 

where d denotes the number of spatial dimensions. The extra spatial coordinate r runs 
from the horizon at r = to the boundary at r = oo, and the temperature of the black 
brane is obtained by demanding the absence of a conical singularity at r^, leading to 
[16] 

T = d± ^ 1 (r h ) z ■ (2-2) 

The metric (2.1) enjoys the Lifshitz isometry, which will be inherited by the boundary 
theory as we will show shortly, and which reads 

r -»• Ar , t -»• \~ z t , x -»• A _1 x , T -»• X Z T . (2.3) 

Next, we introduce a Dirac fermion \I/ into the above gravitational background. Fol- 
lowing [11, 14], we include a boundary action on an ultra-violet (UV) cut-off surface at 
r = vq. While many of the results can be phrased for arbitrary dimensions, for later 
purposes we specify to d = 4 spatial dimensions already here. In the conventions of [14], 
the total action for the bulk fermion that we consider is 



- f d 4 xV^hh{Zlp z (r,x)^ + +ig f ^^\ , 



(2.4) 



where the first term on the right-hand side is the Dirac action in the bulk and the second 
term is the boundary action. Our notation is such that ^ = ^T- with V- being anti- 
hermitian, and a natural choice for the Dirac matrices employed in this paper is given in 
appendix A. 2. Furthermore, we have decomposed ^ into chiral components according 
to its eigenvalue under V- [8, 9] in the following way 
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The normalization constants gj and Z appearing in (2.4) are left unspecified for the 
moment. The coupling of the fermions to the gravitational field is through the viel- 
beins ej 1 and the spin connection (f^)a&, which are determined from the geometry and 
whose explicit form can be found in appendix A. 2. The symbol l}) z appearing in the 
boundary action is the usual Dirac operator for a fermion in a Lifshitz background with 
arbitrary dynamic exponent z, and h denotes the determinant of the induced metric on 
the boundary. 

Following [24] , we also define Fourier-transformed spinors on each constant r slice as 
*±(r,x) = J ^^±(r,p)e ! ^' , p, = (-u,k) , (2.6) 

where ui and k denote the frequency and momenta of the plane wave. The Dirac equation 
{Tp — M)* = resulting from (2.4) (see [14] for more details) and in-falling boundary 
conditions that we shall employ at the horizon imply a relation between ip + and ip—. It 
can be expressed as 

i/j_(r,p) = -i£(r,p)ip + (r,p) , (2.7) 

and can be used to integrate out ip- from the action. Together with (2.4) and (2.6) this 
results in a holographic effective action for the field ip + on the cut-off surface as 

S eS [V + ] = -[_ r [Z0 z (p) + 9fVF"ttro, P )] </>+ • (2-8) 

The Green's function of ip+ that follows from (2.8) for our geometry (2.1) reads 

G R (r ,p) = - (r z V(r )lp z (p) + 9 -l r l + z V 2 (r ) £(ro,p)) _1 , (2.9) 

where we have rescaled ip + so that it acquires a canonically normalized kinetic term. 
The final step is now to take a double-scaling limit 

r -»• oo , g f -»• , 9f r l +Z ~ 2M = const. (2.10) 

As we take this limit, we make sure that the resulting effective action has a kinetic term 
expected from a theory with dynamical scaling z. In the case z 7^ 1, this can be achieved 
by renormalizing away a relativistic kinetic term and adding appropriate counter-terms 
to the action [14]. The result for z / 1 then is 

G R (k, u) = -(u>- jff-k l^r 1 - E(p)) , (2.11) 

where A is an arbitrary real number in our approach, 3 and where in the following we 
denote k = \k\. The above formula for the retarded Green's function is also valid for 

3 In the notation of [14], the parameter A appears as A = — 1/77. 
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z = 1, but then the relativistic invariance requires that A = 1. The self-energy is 
the interaction term that arises from coupling the elementary field ip+ to the conformal 
field theory encoded in the Lifshitz background. As can be seen by comparing with (2.9), 
it is related to the quantity £(r,p) introduced in (2.7) by 

£(p) = -g lim r 2 M Z(r ,p) , —\<M<\, < 2 .12) 

7*0— >-00 v ' 

where g is the coupling constant which stays finite in the limit (2.10) and reads 

9= | r l +Z -2M_ (213) 

The definition of the self-energy in (2.12) is valid for all values of the momentum k, 
however, for k = the allowed range of M is extended to — z/2 < M < z/2. This can be 
derived from the asymptotic behavior of £(r,p) near the boundary, and for more details 
we refer the reader to appendix B. Furthermore, it is possible to extend the relation 
between £(r,p) and even to \M\ > 1/2 for non-zero momentum, which requires 
introducing certain counter-terms on the cut-off surface at ro before taking the limit 
ro — > oo. Since we will also be interested in the range 1/2 < M < z/2, using the results 
in appendix B, we can show that the relation (2.12) in this case should be modified to 

E(p) = - 5r lim (V^cP) - rg^" 1 ) , \ < M < f . (2.14) 

Note that the second term in this expression, which is divergent for M > 1/2, removes 
the divergence in £(r,p) and yields a finite result for 

The transfer matrix £(r,p) defined in equation (2.7) is a complex two-by-two matrix. 
In the case we are interested in, that is, d = 4, it can be diagonalized by choosing 
the Weyl basis of gamma matrices (A. 15). A first-order differential equation for the 
eigenvalues £±(r,p) of £(r,p) can be derived, which was achieved in the anti-de-Sitter 
case {z = 1) in [24] and generalized to arbitrary z in [14]. The resulting differential 
equations read 

r 2 WU ± + 2Mre ± = ^ T h + ± h) & , (2.15) 

where we used rotational invariance to set k = (0,0, £3). Imposing in-falling boundary 
conditions at the black-brane horizon, which corresponds to considering the retarded 
Green's function, leads to the boundary condition 

e±(r h ) = i. (2.16) 

The functions £±(r,p) as well as the resulting self-energy S(p) have various symmetries, 
which are discussed in more detail in appendix C. 
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Figure 2: Illustration of our holographic construction of the self-energy. 



From the imaginary part of the retarded Green's function Gr, the total spectral 
weight, or spectral function, can be obtained, which is of great importance in condensed- 
matter physics and is also directly observable in experiments. It is defined as the trace 
of the two-by-two matrix Gr, i.e., 



It was shown in [14] that the spectral-weight function corresponding to (2.11) satisfies 
the sum rule 



in the physically allowed range for the scale dimensions of the CFT operator coupled to 
ip. In [14], this range was found to be — z/2 < M < z/2, where the lower bound is due 
to unitarity in the field theory, and the upper bound is a consequence of the requirement 
that the coupling of the elementary fermion to the CFT is irrelevant in the ultra-violet. 
Within this allowed physical range for M, the Kramers-Kronig relations are satisfied 
and the sum rule (2.18) is obeyed. The latter is a consequence of the fact that Gr 
is the Green's function for an elementary field ip and its hermitian conjugate that 
satisfy anti-commutation relations. In conclusion, the prescription (2.11) thus allows us 
to holographically compute single-particle Green's functions, which are important from 
a condensed-matter perspective. 

To close this section, we summarize and illustrate our holographic construction in 
figure 2. In particular, equation (2.11) describes chiral single-particle propagators on the 
boundary of the spacetime. These chiral fermions interact with a conformal field theory 
due to a fermion of opposite chirality that travels into the bulk spacetime, has classical 
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interactions there, and comes back to the boundary. In the ultra-violet, the points where 
the chiral fermion emits and reabsorbs the fermion of opposite chirality are very close 
together and the latter fermion cannot travel far into the bulk. Hence, it only notices 
the flat background on the cut-off surface in the bulk spacetime, which leads to a free 
fermionic propagator in the UV. For the infra-red (IR) case, however, the fermion of 
opposite chirality can indeed travel far into the curved part of the bulk spacetime, and 
as a result the IR dynamics will be dominated by interactions and will generally not be 
free. This corresponds physically to the usual "energy-scale" interpretation of the extra 
spatial direction. 



3 Single-particle spectra 

In the previous section, we have outlined the construction of the retarded single-particle 
Green's function. In the present section, we now study this Green's function for both 
the relativistic z = 1 and non-relativistic z = 2 case in more detail, and work out its 
physical properties. The zero-temperature results were first obtained in [14], and here 
we consider non-zero temperatures as well. However, let us stress that the systems we 
describe are always at zero chemical potential. 



3.1 Relativistic case z — 1 

In the case that the elementary fermion ifj + interacts with a relativistic CFT, the back- 
ground is given by an AdS black brane and is described by the metric shown in equation 
(2.1) for z = 1. The mass range for the bulk Dirac fermion is then restricted to lie within 
the interval -1/2 < M < 1/2. 



3.1.1 Zero temperature 

For vanishing temperature and z = 1, the self-energy T,(p) can be computed analytically 
for arbitrary frequencies uj and momenta k. The result for the full retarded Green's 
function has been worked out in [14] and reads 

G R (k,u) = y- — r— -^(u + ff-k), (3.1) 

p2( 1 _ 5Cie -«r(M+i) p 2M-lj V ) 

where we defined p = \/ ui 2 — k ■ k as well as the constant (for arbitrary values of z) 

P (1 M\ 

Note that (3.1) can be extended into the complex plane by allowing for complex momenta 
p, which is important when determining the pole structure of the Green's function. 
Because of the branch point at p = 0, we need to introduce a branch cut which is taken 
to run from p = to p = —ioo for later convenience. Furthermore, we find that in 
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order for (3.1) to be free of singularities in the upper half uj plane, i.e., to satisfy the 
Kramers-Kronig relations, we have to demand 

9 > , (3.3) 

which is derived in detail in appendix D.l. It is then straightforward to show that 
the Green's function (3.1) also satisfies the sum rule (2.18), as it was done in [14] and 
summarized in appendices D.2 and D.3. Furthermore, as mentioned above, (3.1) is valid 
on the complex p plane, with the prescription that on the real uj line the self-energy is 
found by using 

Voj 2 — k ■ k , uj > +\k\ , 

V oj 2 — k ■ k , uj < — \k\ , 

\Jk -k - uj 2 , -|fc| < uj < +\k\ . 

(3.4) 

This directly follows from (3.1) by noting that the region uj < — \k\ for real frequencies 
is obtained from the region uj > \k\ by p — > e lw p, that is uj — > e l7T uj and k — > e ln k, while 
the region \uj\ < \k\ is obtained via p — > e ln / 2 p. Note also that the Green's function for 
— \k\ < uj < +\k\ has no imaginary part. 

3.1.2 Non-zero temperature 

We have studied the effects of non-zero temperature on the structure and form of the 
Green's function by numerical integration of the Dirac equation. To illustrate our results, 
we have included figures 3 and 4 which show the total spectral- weight function (2.17) for 
T = 1/30 and T = 2, respectively. Furthermore, we made a distinction between positive 
and negative masses M, and we employed the rotational symmetry to set k = (0, 0, k%). 
Let us discuss these figures in some more detail: 

• First, we note that a non- vanishing temperature results in a smearing out of the 
features of the Green's function, as can be seen by comparing figures 3 and 4. 

• For M > we observe a large spectral weight approximately at uj = ±£3, which 
is due to the "light cone" of relativistic physics. In the zero-temperature case 
(3.4), the imaginary part of the Green's function is strictly zero outside of the 
light cone in the region — |fc| < uj < +\k\. For non-zero temperatures, there is a 
small contribution in this region, which increases with temperature. 

Note that the peak in the spectral weight is not caused by a pole in the Green's 
function, but by the aforementioned branch cut. Therefore, it has no interpretation 
as a well-defined quasi-particle excitation, but it is still possible to approximately 
determine the dispersion relation, as it is shown in the next subsection. 

• For negative masses M < 0, the spectrum has a similar form. Indeed, there is 
again a large peak at uj = ±£3. However, for low T, the imaginary part of the 



p 2M-l e -m{M+\) 



_ p 2M-l e -i-K(M+\) 
_p2M-l e +ra-(M+±) 



-p 



2M-1 



p 
P 
P 
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(a) M = + 1/4 (b) M = -1/4 

Figure 3: Total spectral function for T = 1/30, g = 1 and M = ±1/4. 



(a) M = +1/4 (b) M = -1/4 

Figure 4: Total spectral function for T = 2, g = 1 and M = ±1/4. 
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Green's function goes to zero at k% = as u — > 0. It has a maximum in between, 
which leads to a broad maximum in the spectrum. Its approximate location is 
determined analytically in section 3.1.3. 

We have also investigated the ratio of the location of the maximum and its width 
as a function of g. This ratio is approximately constant, which means that its 
shape becomes broader as we increase its position, and vice versa. The spectral 
function at = looks very similar to what is plotted in figure 7 for the non- 
relativistic case with z = 2. It is not possible to tune the parameters such that the 
peak is located sufficiently far away from the origin at oj = while its maximum 
remains sharp. Therefore, we cannot interpret this feature in the spectral function 
as a massive particle. Again, we give some further analytic arguments of this 
statement at zero temperature in section 3.1.3. 



3.1.3 Dispersion relation 



In the last subsection 3.1.2, we have illustrated that the spectral weight is peaked at 
particular functions u>(k). In this section, we now investigate the corresponding disper- 
sion relation of the theory. For vanishing temperature this can be done analytically via 
the single-particle Green's function by solving 



ReG^-(k,uj) = 



(3.5) 



By definition, the solutions to this equation determine the dispersion relation of the 
would-be (quasi-)particle. The interpretation as a particle only becomes justified if the 
width is small compared to the energy of the particle. Indeed, a large peak in the 



spectral density function is only obtained when both Re(G R 
as follows from the identity 



(ReG-^ + tlmE) 2 ' 

Let us first solve (3.5). Using the explicit expression (3.1) 
gives two possible dispersion relations, namely 



and (ImX) are minimal, 
(3.6) 

equation (3.5) generically 



id 



±\k\ 



and 



g ciRe 



e -m(M+iyM-l 



1 



(3.7) 



The first solution describes a free, massless relativistic excitation which is always present. 
The second equation in (3.7) should be studied separately in the two regions \uj\ < \k\ 
and \uj\ > \k\. In the first case, using (3.4) we obtain the solution 



P 



1 

gci 



2 

2M-1 



(3.8) 



Since gc\ > and the left-hand side of (3.8) should be real, we conclude that there is 
no solution for \oj\ < \k\. In the second case, \oj\ > \k\, the solution reads 



(3.9) 
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For < M < 1/2, we do not find a solution since gc\ > 0. On the other hand, for 
-1/2 < M < 0, we obtain 



V 



uj = ±Y k ■ k + m 2 , 



m 2 = (gcicos [vr(M + ±)]) 



2 

1-2M 



(3.10) 



which is indeed close to the locations of the maxima in the spectral- weight function that 
were found for non-zero temperature. For zero temperature, equation (3.10) yields the 
exact result. 

Now we can look at the width, by computing the imaginary part of the self-energy 
at the values (3.10). A straightforward calculation shows that, at zero temperature, we 
have 

ImS = -(w-a-fc)tan tt(m+^\ , (3.11) 

with co given by (3.10). However, for momenta small compared to the gap m (in the 
restframe of the would-be particle), we can approximate 



Im £ ~ m tan 



ir Ml 



1 



(3.12) 



The width is then comparable to the gap, and therefore the peak in the spectral density 
does not have a quasi-particle interpretation. One might try to make the width smaller 
by taking the value of M close to the unitarity bound, M = —1/2. However, notice 
that then the gap also narrows down. We conclude therefore that no true quasi-particles 
exist. 



3.2 Non-relativistic case z / 1 

In principle, every value of the dynamical exponent z can be considered using the pre- 
scription shown in (2.11). However, for our purposes we are mostly interested in z = 2 
where the elementary fermion interacts with a CFT exhibiting z = 2 Lifshitz scale 
invar iance. 



3.2.1 Zero temperature 

Unlike the relativistic case studied above, for an arbitrary dynamical exponent z we were 

— * — * 

not able to obtain Gr(oj, k) analytically for the both ui and k non- vanishing. But it is 
possible to determine the zero-temperature result for k = and uj = separately, and 
then restrict the general expression to a great extent. 

Let us therefore start with the case k = and T = 0. From [14] we recall the 
expression for the Green's function as 

G r (6,uj) = m , (3.13) 

uj - gc z uj z e { z + 2> 
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where we employed the definition (3.2). In the case with to = and T = we instead 
find [14] 

0) = i - r, ? i 1 ,2M t~ T ' ( 3 - 14 ) 

f a ■ k k z ~ l +gci k 2M ~ 1 a ■ k v ; 

where again k denotes \k\. The generic case at vanishing temperature when both k 
and oj are non-zero can then be restricted as follows. Using the scaling and rotational 
symmetries of the self-energy (see appendix C.l), and defining u = ui/k z for notational 
simplicity, we have 

Gn(k,u}) = ^- , (Ql^ 

u - \a ■ k k z ~ l - g (k 2M s 1)M («) + k 2M - 1 s 2 , M («) a ■ kj { ' 

where s± t M and «2,M are complex functions of u and M only. Furthermore, we can 
derive conditions on the functions s\p using the symmetry of S(p) under the change 
M —7- —M, as discussed in appendix C.4. The asymptotic behavior of the functions 
si,m(u) and S2,m(u) in the limits u — > and u — > oo is fixed by the expressions (3.13) 
and (3.14) given above. We therefore have 

si,m{u) > u*e \z^2) Cz , s 2 ,m{u) > 0, . . 

/ \ u^O n , s u^O {6.10) 

si,m{u) > , s 2 ,m{u) > Ci . 

Unfortunately, these conditions do not seem to be sufficient to determine the analytic 
form of (3.15). Therefore we have to study the Green's function numerically. Yet, we 
can determine the qualitative form of the dispersion relation by analytic arguments as 
we describe below in section 3.2.3. 



3.2.2 Non-zero temperature 

We have studied the retarded Green's function and the corresponding spectral-weight 
function for z = 2 numerically as a function of w, k, T and A. Using the symmetries 
summarized in appendix C, that is chirality and particle-hole symmetry, we observe that 
the components of the spectral-weight function obey the following relations 

p ± (fc,w) = p T (-k,u) = p ± (-k,-uj) . (3.17) 

Consequently, it suffices to consider separate components instead of the trace over chiral 
components. 

Next, we recall that for z ^ 1 the Green's function (2.11) contains the parameter A, 
which we have not determined analytically. Treating A as a free parameter, we observe 
that under a change of sign of A, the far UV behavior of one component of the spectral 
weight function asymptotes the UV behavior of the other component with the original 
sign, that is, 

Pt x &») M " 00 ' |fcK0 °: p% x & U ) = (3-18) 
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Since A can be interpreted as a spin-orbit coupling constant, we use the following con- 
vention for plotting components of the spectral-weight function: when A is positive 
(negative), the plus(minus) component of the spectral- weight function is shown. In this 
manner we make sure that we always compare spectra with equal group velocity in the 
UV, also when A changes sign. We make this choice because qualitatively the UV physics 
is then, apart from the topology of the band structure, independent of the sign of A. 
This allows us to compare more clearly the physics for positive and negative values of A 
as we will see shortly. 

In figures 5 and 6 we show numerical results for separate components of the corre- 
sponding spectral-weight function. The parameters are chosen as M = ±1/4, g = 1 and 
T = 1/30, for a number of different values for A. We discuss some of the features in 
turn: 

• In figure 5, taken for positive M, we see that the spectral-weight function is very 
sharply peaked in the UV and behaves as a free chiral fermion with a quadratic 
dispersion. In the IR, the self-energy becomes dominant which changes and smears 
out the form of the free-fermion dispersion. When A switches sign, the spectral- 
weight function changes in the IR. However, note that our convention is such that 
the UV behavior stays the same. 

• In figure 6, where M is negative, we observe a phenomenon similar to the relativistic 
case with M < 0. In particular, the spectral weight has to vanish for k = 
as co — > 0, so there is a maximum in the spectral-weight function at non-zero u. 
However, the ratio of the location and the width of the peak again remains constant 
as we change g, so there is no gap generation with a quasi-particle interpretation. 
This is further illustrated in figure 7. 

To close our discussion about non-zero temperature effects in the case of z = 2, let us 
consider the strict IR (or hydrodynamic) limit u) — > 0, k — > 0. As explained in appendix 
B.3, in this limit the contribution from the free propagator vanishes and the Green's 
function reduces to the inverse of the self-energy An analytical result for all values 

of the dynamical exponent z and dimension d and for non-zero temperature has been 
obtained for this case in equation (B.21), which we recall here for convenience 



3.2.3 Dispersion relation 

We now consider the dispersion relation for z ^ 1 at zero temperature, which can be 
derived from the general form of the Green's function given in (3.15). In particular, the 
dispersion relation is obtained by solving (3.5), which in the present case reads 



2M 




(3.19) 



u - - B ■ k k z - 1 - g k 2M Re s lM («) + 



a ■ k 



S2,M (u) = . 



(3.20) 



k 
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(e) Density plot of p (k, ui) at A = —0.5 (f ) Density plot of p (k, to) at A = 

Figure 5: Density plot of p ± (k, oj) for z = 2, M = +1/4, T = 1/30, g = 1. 
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(e) Density plot of p (k, ui) at A = —0.5 (f ) Density plot of p (k, to) at A = 

Figure 6: Density plot of p ± (k, oj) for z = 2, M = -1/4, T = 1/30, g = 1. 
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p ± {0,uj) k = 0, T = 0, M = -1/4, <q = 0.4574 




Figure 7: Spectral- weight function for vanishing momentum determined from (3.13). Since the 
ratio of the location and the width of the peak remains constant, there is no quasi-particle 
interpretation. 



In the following we consider the upper-spin component for definiteness. The latter means 
that a ■ k/k is replaced by sign^), where we again employed the rotational symmetry 
to align the momenta in the z-direction. For positive (negative) momentum k% the 
dispersion relation then becomes 



1 iz 



g k 2M Re (s 1|M («) ± s 2 ,m («)) = . 



(3.21) 



In order to determine the qualitative form of the dispersion relation, we study various 
limits of the equation above. 

• We first consider the UV limit u — > 00, k — > 00 for which there are the three 
distinct possibilities: uj/k z — > 00, k z /ca — > 00 and oj/k z — > const. The first two do 
not allow for a solution of (3.21), but the third possibility leads to 



4- 



for k > 1 



(3.22) 



in the allowed mass range — z/2 < M < z/2. This result confirms our general 
picture that the interaction of the elementary fermion with the CFT is irrelevant 
in the UV, hence the dispersion becomes that of a free fermion. 

Next, we consider the limit k% — > + , oj > 0. In the mass range —z/2 < M < z/2 
we determine the dispersion curve from (3.13) as 



u + gc z 



cos 



sign(M) u>~ 



. 



(3.23) 
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Noting that g > and c z > 0, we find that the only solution to this equation is 
oj = + for M > 0. On the other hand, for M < we obtain the only solution as 



These solutions correspond to the points where the dispersion curve touches the 
axis at k = 0. 

Studying the dispersion relation for the case of ks — > + and oj < 0, we find 
another solution at w = —m z , consistent with the particle- hole symmetry. This 
is the non-relativistic analog of the result we obtained in equation (3.10). As 
mentioned above, and similar to the case with z = 1, the corresponding peaks in 
the spectral-weight function cannot be interpreted as a true gapped quasi-particle 
excitation. This is further illustrated by a plot of the analytic result in figure 7. 

• The most interesting case is when < M < 1/2 and A < 0. As one can see 
from figure 5, this situation is quite different from < M < 1/2 and A > 0. The 
dispersion curve oj(k) for A < crosses the oj = axis three times, instead of 
only once at oj = k = 0. It therefore seems that two Fermi surfaces appear at the 
momenta for which oj(±kp) = 0. This signals a phase transition as one crosses 
from positive to negative A. The remainder of this paper is devoted to a detailed 
study of this phase transition. 

4 Quantum phase transition 

In figures 5 and 6 we have already seen that for z = 2 the parameter A controls features of 
the spectral density significantly. In particular, when switching the sign of A the number 
of zeros of the dispersion relation oj(k) changes, as becomes clear when comparing for 
instance figures 5(a) and 5(f). In the present section, we study this phenomenon in 
more detail. For simplicity we focus on the case < M < 1/2 since the quantum phase 
transition discussed here is also present for —1/2 < M < 0. However, we also briefly 
consider the case \M\ > 1/2 in section 4.5. 

4.1 Momentum distribution 

Let us start by defining the momentum distribution function. We choose again a con- 
vention such that the group velocity of a particle with one of the spin components always 
has the same sign at large momenta, irrespectively of the sign of A. More concretely, 
with uf(oj) the Fermi distribution we write 




M < . 



(3.24) 




for A > , 



for A < . 



(4.1) 
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The behavior of N- when A changes sign is illustrated in figure 8 for non-zero tem- 

k _ 

peratures, where we again employed the rotational symmetry to set k = (0,0, £13). In 
particular, for A < the momentum distribution indicates a Fermi surface with a cer- 
tain width, on which we comment later. When A switches sign, two extrema appear that 
develop into sharp discontinuities as A increases. This is a clear signature of a Fermi 
surface. The locations of the jumps are determined analytically in the next section. 

4.2 Fermi momentum 

We now study analytically how the number of putative Fermi surfaces changes as we 
vary the parameter A. To do so, we recall that Fermi surfaces are determined by the 
poles of the Green's function at vanishing frequencies oj = 0. The latter can be computed 
by setting to zero the denominator in (2.11). Note also that the self-energy is given by 
(2.12) which satisfies the particle-hole symmetry derived from (C.16). Employing the 
general form of the self-energy implied by (3.15), we see that in the limit oj — > the 
imaginary part of S(p) vanishes and the condition for the presence of a Fermi surface 
coincides with the presence of zeros of the dispersion relation (3.5). Thus, the loci of the 
Fermi surfaces are given precisely by the points where the dispersion curve oj(k) crosses 
the oj = axis. 

We therefore consider equation (3.21) in the limit oj — > 0. Using (3.16) and the 
expression (3.2) for ci, we obtain the following formula for the loci of the Fermi surfaces 

\k z + g k 2M Cl = . (4.2) 
A 

Then, as we have mentioned in equation (3.3), to avoid violation of causality we have to 
require g > 0. Equation (4.2) can thus have non-trivial solutions in the range —1/2 < 
M < 1/2 only when A < 0, that is 

k = k F ={-g\c l )^ B . (4.3) 

Next, let us consider again the numerical results for the imaginary part of the Green's 
function shown in figure 5. For A < 0, there are zero-energy modes at non-vanishing 
momentum, suggesting that the system indeed has two Fermi surfaces at zero temper- 
ature. However, due to the small but non-zero temperature which we have to employ 
in our numerics, the spectrum shown in the plots has a finite width and the locations 
of the zero modes are approximately at the Fermi momentum (4.3). To investigate this 
point further and to confirm that we are indeed dealing with a genuine Fermi surface, in 
the following we compute the quasi-particle weight Z as a function of A and g, as well 
as the effective mass and the lifetime at the Fermi surface as a function of oj, T, g and 
A. 

4.2.1 Quasi-particle residue 

In order to scrutinize the Fermi surfaces, in this subsection we determine the quasi- 
particle residue by linearizing the dispersion relation around the Fermi surface at oj ~ 0, 
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(e) A = -1, k F = 1.63 (f) A = -100, k F = 35.2 



Figure 8: Momentum distribution of the minus component for M = +1/4, T = 1/30, g = 1. For 
A > it has a smooth kink-like behavior whereas for A < two Fermi surfaces develop. The 
red dots give the analytic value (4.3) of kp. For small and negative A, the Fermi surfaces are 
smeared out because of the non-zero value of the temperature. 
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k ~ k F . For convenience, we choose the lower component of the Green's function 
which results in the spectral- weight function shown in figures 5(d)-(f). Furthermore, we 
employ the rotational symmetry to set k = (0, 0, k^) with k% > 0, and thus k^ = k in the 
following. We then compute up to first order in derivatives 



Re (G^)- 1 

ww (l-d UJ Re'E-(k F ,uj) J - (jfe - k F ) d k ^ Re 0) - ^k z ^j 



(4.4) 

k=kp 



Using (4.4) in the retarded Green's function, close to the Fermi surface we obtain the 
expression 

—Z 

G R (k,uj) - - _ _ — — ^ ReS _^ )0 ^ _ 1 -iZImE-(k,u) ' ^ 
where the wavefunction renormalization factor, or quasi-particle residue Z, is given by 

Z = l-a w ReS±(fc F ,w)| w=0 ' (4 - 6) 

The residue is equal for both spin components because the ui derivatives of Re S 1 * 1 are 
equal at oj = 0. Furthermore, note that Z depends on the coupling g explicitly due to 
the factor g in S(p), but also implicitly via the dependence of kp on g shown in (4.3). 
The dependence on A is only through k F . To be able to determine the quasi-particle 
residue, the real part of the self-energy has to be linear in oj = 0, which is indeed the 
shown in figure 9. The first derivative of the real part of the self-energy is shown 
for T = 1/100 in figure 10(a). 

Let us illustrate the calculation of Z for a particular value of A and g, namely 
A = —0.5 and g = 1. For T = 1/30, the first derivative of the real part of the self- 
energy has the non-zero value <9 W Re T* + (kp, 0) « —0.8 at u = 0. This leads to a finite 
quasi-particle residue of about Z 0.56, which is precisely the height of the step in the 
momentum distribution at ks = ±fcp, as shown in figure 11. Repeating this calculation 
for different values of A, we have determined Z{\), which is shown in figure 12. Note 
that < Z < 1 as is required. 



4.2.2 Effective mass 

In the last subsection we have seen that at the Fermi surfaces there are quasi-particles 
which are strongly renormalized. Now, to further characterize the Fermi liquid we should 
also compute their effective mass. The latter is defined as the inverse of the slope of the 
quasi-particle dispersion around the Fermi surface, that is, 

uj{k) = (k — k F ) k F near the Fermi surface, (4.7) 
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Figure 9: Imaginary and real part of the self-energy S + (p) evaluated at the Fermi momentum 
for T = 1/100, g = 1, A = —0.5 and M = 1/4. The imaginary part is zero around lu = and the 
real part is linear around ui = 0, which are both defining properties of a Fermi liquid. 



a. Re 




-1.2 L 

(a) cL Re E 4 



(b) dllm^ 



Figure 10: This figure shows the first derivative of Re£ + (p) and the second derivative of 
Im£+(p), both at A = -0.5, T = 1/100, M = 1/4 and .g = 1, evaluated at the Fermi mo- 
mentum. 
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A = -0.5 
k F = 1.0305 
Z = 0.557 




Figure 11: The momentum distribution for M — 1/4, T = 1/30 and g = 1. The quasi-particle 
residue is given by the distance between the dotted lines: the difference in height is the numerical 
value of Z, which is 0.56 for A = —0.5. The red dots give the analytic value of ±fcp-. 
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Figure 12: Wavefunction renormalization factor Z as a function of A for g = 1, M = 1/4 and for 
different temperatures. Z is only defined for negative A. The quasi-particle weight lies between 
one and zero, as expected. It approaches unity for very large and negative values of A, i.e., very 
far away from the phase transition. For A — ^ — it vanishes only strictly at T = 0. Due to 
non-zero temperature effects, the quasi-particle residue has a non-zero minimum at A = 0, which 
is smaller for the lower temperature curve. The inset shows a plot of m c ff/|A| as a function of A 
for M = 1/4. 
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where we are again employing k^ > and thus ks = k. Therefore, in our present 
situation the retarded Green's function near the Fermi momentum can be written as 

— Z 

Gr{K W) = cj-^(k-k F )k F -iZJmE-(k,u) ' (48) 
and by comparing with the explicit form given in (4.5) we conclude that 



— = ^8 k (ReE-(M)-T^ 
m cS k F V A 



k=kp 



(4.9) 



At zero temperature, we can calculate the effective mass analytically. In particular, for 
z = 2we obtain from the exact result given in equation (3.14) that 

ReZ-(k,0)=gc l k 2M , d k Re S"(fc, 0)\ k=kp = ~ k F , (4.10) 

where we eliminated gc\ in favor of k F . Therefore, the effective mass for z = 2 at zero 
temperature reads 

meS = -2Z(X)(M + l)- (4 " U) 
The effective mass is shown in the inset in figure 12. 

4.2.3 Quasi-particle decay rate 

We now turn to the lifetime of the quasi-particles. At zero temperature, from the 
analytic result (3.14) we infer that the imaginary part of the self-energy vanishes at 
u = 0, k = k F . In figure 9, which shows the imaginary part of the self-energy, we see 
that this behavior is confirmed numerically. This implies that the spectral function at 
the Fermi surface is a 5-function at T = 0, which is a well-known and defining property 
of a Fermi liquid. As a consequence, the momentum distribution shows a discontinuity 
at the Fermi surfaces. Now, the inverse of the imaginary part of the self-energy Im S 1 * 1 
has SI units of time, and at co = 0, we associate it to a lifetime for the quasi-particles 4 

1 

T k = T7^ : ■ (4.12) 

fc -2ZIm£±(£;,0) 

The lifetime T kp is infinite precisely at the Fermi momentum. We can also define a 
frequency-dependent decay rate at the Fermi momentum as 

r(w) = -2ZImS ± (fc F ,w) . (4.13) 

At the Fermi surface oj = 0, the decay rate T is the inverse of the quasi-particle lifetime, 
therefore it vanishes for oj — > 0. This behavior can also be observed in the spectra above. 



4 Here, we ignore cj-dependent corrections to the lifetime, as the dominant contribution comes from 
the behavior at lo — 0. 
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g = 0.25 blue curve 
g = 0.5 red curve 
g = 1 yellow curve 




(a) T= 1/20 (b) T= 1/50 



Figure 13: The spectral function for the spin-up component at the Fermi momentum as a function 
of u. In both plots parameters are chosen as gX — —0.5. This figure confirms numerically two 
important points that were made in the text. Firstly, for the higher temperature (a), a sharp 
peak is just visible at u = 0. For the lower temperature (b), this sharp peak is still there but is so 
narrow that it can barely be discerned. This indeed indicates that, at T = 0, there is a (^-function 
behavior at the Fermi surface, and that this (S-function peak becomes broader as temperature 
increases. Secondly, close to the Fermi surface, that is, for small Li, the form of the spectrum 
depends only on the product gX. Indeed, all four curves for different g but equal gX overlap for 
small ui. This overlap is better for lower T. This shows that the scale of the Fermi momentum 
is indeed set by the combination gX and that the spectra are close to the Fermi surface indeed 
governed by kp from (4.3) as we expect. 



In particular, when T and lj are small but non-zero, the <5-function is broadened, which 
is indeed visible numerically, as shown in figure 13. For completeness, we note that the 
retarded Green's function given in (4.5) can now be written as 

^ = — i ft.Zl.t.^i ■ ( 4 - 14 ) 



LJ 



(k - k F )k F + ^ 



A broadening of the (^-function implies a non-zero decay rate for the quasi-particles. 
For a normal Fermi liquid, the imaginary part of the self-energy at the Fermi momen- 
tum vanishes at zero temperature and small frequency as lj 2 , which stems from Pauli 
blocking for scattering near the Fermi surface. For small non-zero temperatures, thermal 
averaging over lj 2 yields a T 2 dependence of the decay rate. 

Let us compare this with our numerical results. We observe that the decay rate 
vanishes much faster than in a conventional Fermi liquid. In figure 10, we show the 
first derivative of Re S + (p) and the second derivative of Im E + (p) as a function of w 
at k = kp. The former is non-zero at lj = 0, which indicates a non-zero value of Z as 
explained above. The latter is zero at lj = 0, which shows that the quasi-particle decay 
rate vanishes faster than to 2 . The behavior shown in figure 10 is in agreement with the 
existence of Fermi surfaces at k = kp, however the decay rate is more suppressed than the 
typical oj 2 behavior of a conventional Fermi liquid. This indicates that, next to Pauli 
blocking, there is another mechanism at hand which strongly suppresses interactions 



27 



around the Fermi surface for low temperatures T and small frequencies oj. This appears 
to be a general feature of holographic Fermi liquids and can be understood better by 
performing a WKB analysis of the bulk Dirac equation. 

4.3 Second-order WKB analysis around the Fermi surface 

In section 4.2.3 we have seen that the fermionic quasi-particles are rather stable since 
their lifetime at the Fermi surface diverges faster than the power law expected from 
Fermi-liquid theory. To gain more insight into this property, we would like to obtain an 
analytic expression for the self-energy at the Fermi surface, that is, we want to know 
the form of the self-energy for both oj 7^ and k = kp ^ 0. As stated before, it is 
not possible to solve the differential equation for £ exactly for generic non-zero values 
of oj and k, however, we can deploy the WKB approximation to study the fermionic 
excitations around the Fermi surfaces analytically, as was done for instance in [25, 26]. 

In principle, the WKB calculation gives the exact result for the retarded Green's 
function in the limit k — > kp and oj <C 1. Indeed, as we will see, it yields very elegantly 
the functional form of the self-energy at the Fermi surface, from which we can derive 
the behavior of the quasi-particle lifetime at k = kp. The behavior turns out to be 
exponential, i.e., 

ImS~ e - a ^, (4.15) 

with a a constant to be determined. This result is consistent with the exponential be- 
havior found in [25, 26]. Furthermore, we also find, to first approximation, an analytical 
result for the wavefunction renormalization Z, which is consistent with the strong quasi- 
particle renormalization found numerically. Notably, to obtain the latter result, we have 
to use the second-order WKB approximation. 

However, before presenting the explicit result for the self-energy of our second-order 
WKB calculation, we have to mention the following. Obtaining the exact prefactors 
of the various terms in the final result, e.g. the constant a in (4.15), turns out to be 
non-trivial. The reason is that several integrals in the computation cannot be performed 
analytically (see appendix E for details) but have to be evaluated numerically. Alter- 
natively, we could make an approximation to the integrals which, however, turns out 
to be similarly difficult. Our approach here is to deploy the WKB approximation to 
gain an understanding of the functional form of the quasi-particle lifetime at the Fermi 
surfaces. Therefore, the result that we present in equation (4.21) is obtained using the 
latter method, and so the prefactors, including a, are not exact. At the end of appendix 
E, we show some numerical evidence for the conjectured value of the constant a. 

Let us outline the basic steps of our WKB calculation and refer to appendix E for 
further details. We determine the Green's function at zero temperature. The expansion 
parameter that plays the role of 1/h in the conventional WKB approximation, and that 
we take to be large in order to satisfy the WKB conditions, is the rescaled momentum 
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(see appendix C.l) 5 



k 

> 1 . (4.16) 



To simplify the calculation we also assume that both u> and k are positive, and we 
introduce the rescaled variable x = r/\fuj. Choosing the up and down components u± 
and d±, respectively, of ij)± introduced in (2.6) as 

fa = T^W) Ct) ' (4 ' 17) 

we define new fluctuation fields y± in terms of the original fields u± as 



^F l u + . (4-18) 



y± = x s/2 u ± 

Then, the once-iterated Dirac equation can be written as a Schrodinger-like equation 

' " 2 y±(x)-V ± (x)y ± (x) = 0, (4.19) 



k 2 dx 2 

with the effective Schrodinger potential at T = of the form 

1 1 (|±M)(f±M) 3 1 §±M 
±(X) "^"P^ + + 4( x fe±l) 2 " xfc(xfeil) ' (420) 

The characteristic behavior of the Schrodinger potentials for the chiral components is 
illustrated in figure 14. In particular, V + has a single classical turning point at x tp and 
V- has a pole and a turning point. For k S> 1, both effective potentials have a very 
high potential barrier close to x — > 0. For simplicity we present here the calculation only 
for the plus-component, the result for the minus-component is easily obtained using the 
symmetry explained in appendix C. Also, we only perform the calculation for the case 
M > so that the single turning point of V+ at xt p divides the entire range into two 
regions, a classically allowed region where V < and a classically disallowed region 
where V > 0. 

The strategy to determine the WKB solution is as follows. We first incorporate the 
in-falling near-horizon boundary condition into the WKB wavefunction in the classically 
allowed region x < x tp . Then, we apply the usual linearized connection formulae at 
x = x tp to continue the solution into the classically disallowed region x > x tp - Next, we 
match this solution onto the exact near-boundary solution to fix the remaining integra- 
tion constants. Finally, from the asymptotic expansion of the solution we read off the 
coefficients A± (see appendix B.2) that determine Gr via £ = iA_/A + . The details of 



5 We will specialize to the case of interest z — 2 in this section, although the calculation can easily be 
generalized to arbitrary z > 1. 
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(a) Linear plot of the effective potential V± 
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(b) Double-logarithmic plot of \V±\ 



Figure 14: Plots of the effective potentials V+ and V- (blue and red line) as function of x = r/^/uj 
at T — for k — 10 and M = 1/4. In the near- horizon region < x < Xt p , both effective 
potentials are negative. The potential V+ has a turning point at Xt p , which can be seen as the 
first dip from the left in (b) . V_ has a pole located approximately at the position of the classical 
turning point of V + . V + is always positive for x > x tp , while V- has a different turning point at 
an intermediate value of x. This point can be seen in the inset in (a) and as the second dip from 
the left in (b). 



this calculation are explained in appendix E, and the final result for the WKB self-energy 
reads 

1.2 

1 — i 7 ~ c~ ilxM c~ 2 
E+(k,u) =g Cl k™ 2± ^± , (4.21) 

where we employed again k = (0, 0, k$) with ^3 > 0, so that the corresponding Green's 
function in the limit k S> 1 reads 

Gi t $, U ) = -L-\v-g Cl k™ ?± ^ . (4.22) 

The constant c\ was given in (3.2), and we defined 

As mentioned before, the result (4.21) comes from performing an expansion within the 
framework of the WKB calculation. In this case, we find a = 2 for the aforementioned 
coefficient in (4.15), which we can improve by taking into account more terms in the 
expansion. In particular, as explained in more detail in appendix E, the correct prefactor 
in the exponent seems to be approximately a ~ ir/2. 

After having obtained our result, let us now study the Green's function (4.22) in more 
detail. As a first consistency check, we see that in the limit cj — > 0, which corresponds to 
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k = k/y/ui — > oo, we find the same locus for the Fermi surface at k = kp as in equation 
(4.3). Furthermore, the decay rate of the quasi-particle excitations around the Fermi 
surface is calculated using (4.12) and employing the imaginary part of the self-energy at 
k = kp, which reads 

ImS+(fe F)W )= 1 g Cl k 2 F M sin (Tr(M-l))^ (^ + — ) e~ 2 # , (4.24) 

where kp is given by (4.3). We note that j± and e± become unity in the strict Fermi- 
surface limit, and the result thus shows an exponentially suppressed decay rate with 
additional ojjkp corrections. This result demonstrates the extremely stable nature of 
the quasi-particle excitations on the Fermi surface and is consistent with previous WKB 
calculations that similarly find exponentially narrow decay rates [25, 26]. A comparison 
with the numerical results is made in figure 15. 

Another analytic result that we extract from (4.21) is the wavefunction renormaliza- 
tion (4.6) on the Fermi surface. This follows from the real part of the self-energy which 
reads 

Re£+(fc,w) =g Cl k 2M ^± , (4.25) 

up to exponentially small terms multiplied by exp(— ak 2 /oj) that will drop out in the 
limit oj — > after differentiation. Combining (4.25) with (4.3) and (4.6), we arrive at an 
analytic result for the wavefunction renormalization factor of the form 

~ 2A , (4-26) 



2A-1 ' 

which for A < satisfies < Z < 1 as required. We note that, to obtain this non- 
trivial result different from one, we have to go beyond the leading order in the WKB 
approximation, since otherwise the derivative of the self-energy in (4.21) would vanish 
due to the exponential suppression. Interestingly, because all factors M in (4.23) are 
multiplied by oj, these drop out in the limit oj — > so that the analytic result (4.26) is 
independent of M. This is in contrast to the numerical results which do depend on the 
value of M. The analytic result (4.26) is therefore expected to be the M = part of 
a larger expression, that one may obtain using a more sophisticated WKB calculation. 
This presumption is supported by figure 16, which shows that both results are consistent 
in the case M = 0. However, the aim of the WKB calculation was merely to support 
our numerical results, and finding perfect agreement is beyond the scope of this work. 



4.4 Phase diagram 

The results obtained in the sections above, in particular the plots in figure 8, can be 
explained by the existence of a quantum phase transition at A = and g > 0. Firstly, we 
can understand why there must be different phases. Consider again the Green's function 
(2.11) and recall that A can change sign whereas g is required to be positive (in the range 
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Dashed red curve: numerics at T = 1/50 
Dashed blue curve: numerics at T = 1/750 
Solid yellow curve: WKB from (4.21) for a — 2 
Solid green curve: WKB from app. E for a = 5/3 
Solid blue curve: WKB from (4.21) for a = tt/2 



Figure 15: Comparison of WKB and numerical results for |ImS + | at k = kp with A = — 1 
yielding kp ps 1.64, and M — 1/4 as a function of u. The dashed curves are numerical results, 
the solid curves are second-order WKB results at T = 0. Note that the WKB result (4.21) for 
a = tt/2 agrees very well with the low temperature numerical result. 
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Dashed red curve: numerics at T = 1/30 for M = 1/4 
Dashed blue curve: numerics at T = 1/30 for M = 
Solid blue curve: WKB at T = from (4.26) 
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Figure 16: Comparison of WKB and numerical results for the quasi-particle residue Z from (4.6) 
as a function of A. The dashed curves are numerical results, the solid curve is the second-order 
WKB result at T = 0. This figure shows that the WKB result for Z resembles the M = result, 
indicating that the analytic expression may be the M = part of a more general result. 
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—1/2 < M < 1/2) in order to satisfy the Kramers-Kronig relation. Furthermore, the 
kinetic term scales as k z , but the self-energy is proportional to k 2M which is a lower 
power of k. As a consequence, there will always be a small but non-zero k interval at 
low values of the momentum in which the self-energy term dominates. When there is 
a difference in the relative sign, there can be a cancellation between the kinetic and 
the interaction term in the Green's function at low momenta. This is what brings the 
Fermi surfaces into existence. When there is no relative sign difference, the self-energy 
contribution adds up to the kinetic term, which results in a different behavior of 

k 

and therefore in a different phase. Thus, the crucial property is the possibility of a 
relative sign change between the kinetic and interaction term which is dominant at low 
momenta. 

We sketch the phase diagram of our system in figure 17, and relate it to our results 
for the momentum distribution at various values of A and T. 

• For A < and temperatures T small compared to the Fermi energy €f = k F /\X\, we 
have a phase with two Fermi surfaces: one for the particles and one for the holes, 
which can be seen as two sharp discontinuities in the momentum distribution at 
ks = ±\kp\. The system is formally only a Fermi liquid at exactly zero temperature. 
Numerically we can obtain only non-zero temperature results, and the effect of 
temperature is to smear out the discontinuities over an interval around ±|fcp|. 
However, when T/cf is small (outside the critical region) we still see the Fermi- 
liquid-like or degenerate behavior. Indeed in figure 8 for small A < the jumps 
are not infinitely sharp but still very much located around a single value of k. 

• When increasing temperature, in the region T ~ tp there is a crossover between 
the Fermi-liquid phase and the quantum critical phase, indicated by the dashed 
lines in figure 17. To calculate the crossover temperature as a function of A we 
realize that the crossover occurs when T « ep, thus at 

T^e F = j^k z F = \X\^M =gc 1 k 2 F M ■ (4.27) 

In the case M = 1/4, z = 2, the crossover is at T ~ (A) 1 / 3 . 

• When T S> ep the Fermi surfaces will be smeared out to an extend that they are 
difficult to identify, as it can be seen in figure 18(a). Then the system is in the 
quantum critical phase. As A — > 0~ , ep decreases, and as a consequence the Fermi 
surfaces develop into smooth bumps for very small A. Numerically, we will always 
observe temperature effects for small values of |A|. 

• At A = and T = 0, there is a quantum phase transition from a phase with two 
Fermi surfaces to a phase with no Fermi surfaces. At A = and T > 0, the system 
is in the quantum critical phase. Interestingly, there is a qualitative difference 
between the quantum critical A < and A > phases, see figure 18(b) for an 
impression of the behavior of N^ 1 for very small A. 
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Figure 17: Sketch of the phase diagram in the case —1/2 < M < 1/2 and g = 1, showing the 
non-Fermi-liquid phase for low T and positive A, the Fermi-liquid phase for low T and negative 
A, and two qualitatively different quantum-critical phases for high temperatures and A positive 
and negative. The system is strictly only a Fermi liquid at T = (thick line), but has Fermi- 
liquid-likc behavior for low temperatures. The dashed curves correspond to the crossovers to the 
quantum-critical regimes. There is a quantum phase transition at A = 0, T = 0. 

• In the case A > 0, the momentum distribution has a kink-like behavior and is 
centered around k = 0. For T <C e_F> the width of the kink is determined by 
kp = (gXc\) 1 ^ z ~ 2M \ whereas for T 3> cf it is determined by T, see figure 18(c). 

4.5 Fermi surfaces and the quantum phase transition for | < \M\ < | 

In the previous sections, we have studied the system in the particular range < M < 1/2 
and we have demonstrated the appearance of Fermi surfaces for A < 0. However, we 
obtain a similar behavior also for positive values of A if we consider 1/2 < M < z/2 and 
-z/2 < M < -1/2. 

Let us demonstrate this for 1/2 < M < z/2. In this case, the relation between the 
self-energy T,(p) and the bulk quantity £(r,p) is given by equation (2.14). However, the 
expression (C.24) for the Green's function is unchanged because this form is determined 
solely by the symmetry properties of £,(r,p), which are also satisfied by the counter-term 
in (2.14). In particular, the ui — > limit of the dispersion relation is again given by 
equation (4.2) which we recall for convenience, 

\ k z + gk 2M a = . (4.28) 
A 

Now, we can find a non-trivial, real solution to this equation in the range 1/2 > M > z/2 
for positive A because the constant c\ given by (3.2) is negative in this range. Note 
that the same conclusion applies for —z/2 < M < —1/2. Thus, we now obtain Fermi 
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(c) A = +1 



Figure 18: Momentum distribution in the quantum critical region for T = 5, g = 1, M = 1/4 
and various values of A. For A < 0, the momentum distribution always shows two extrema at 
non-zero fc 3 , which develop into Fermi surfaces when the temperature is lowered. This is shown 
in (a). The behavior for small negative A is illustrated in figure (b) with kp ~ 0.026. Note that 
because kp is relatively small, we plotted NZ here as a function of k^/kp for clarity. For A > 
the momentum distribution is very broad as can be seen in (c). 



35 



surfaces for positive A in these ranges for M. Furthermore, we also find a quantum phase 
transition at A = below which the aforementioned Fermi surfaces disappear. 

5 Conclusion and discussion 

In our recent work [14] we proposed the construction of a holographic fermionic retarded 
Green's function modified to describe single-particle correlations, using the AdS/CMT 
correspondence. Two aspects of this prescription that make it particularly interesting for 
condensed-matter applications, are the following. Firstly, the retarded Green's function 
satisfies the zeroth frequency sum rule and the Kramers-Kronig relation both at zero 
and also at non-zero temperatures. We have checked that this is the case analytically 
where possible, see appendix D, and otherwise numerically. This shows that the cor- 
relation function is a genuine single-particle correlation function that can be compared 
to experimental ARPES data. Secondly, it can describe quantum critical points with 
relativistic as well as non-relativistic dynamical exponent z > 1, as a consequence of the 
(approximate) Lifshitz isometry of the used bulk spacetime. 

In this paper, we have investigated the physics arising from this modified holographic 
prescription. The resulting spectral- weight function and dispersion relations that we 
presented in section 3, show that we are describing a particle- hole symmetric, chiral 
boundary system that is a strongly interacting Weyl semimetal. Weyl semimetals have 
recently received considerable interest in the condensed-matter community because of 
their fascinating and unusual properties that result from the topological nature of the 
band structure [27, 28, 29, 30, 31, 32]. In order to investigate the band structure, 
the Weyl semimetals are usually treated as free or weakly interacting systems. Using 
our holographic prescription of Weyl semimetals, we are now also able to explore these 
systems when they are strongly interacting or critical. In particular, the inclusion of 
a holographic self-energy leads to the existence of several different phases including a 
Fermi-liquid phase and a non-Fermi-liquid phase, separated by a quantum phase transi- 
tion, as we have shown in section 4. 

There are various interesting directions in which our prescription can be extended. 
Firstly, our current prescription is intended to describe single-particle correlation func- 
tions, which appear often in condensed matter. However, in experiments one also encoun- 
ters two-particle correlation functions. Important examples include the current-current 
correlation function that determines the electrical conductivity, a, and the heat current- 
heat current correlation function which is related to the thermal conductivity k. The 
electrical conductivity has been investigated in a holographic context before, starting 
with [33] , but it would be interesting to also consider this particular function using our 
holographic prescription. In particular, since our present discussion is at zero doping, 
that is, it contains an equal amount of particles and holes, we expect the corresponding 
conductivity in the boundary system to remain finite in the ui — > limit even without 
impurities. Additionally, in Fermi- liquid theory the thermal conductivity n scales as 
1/T, whereas the electrical conductivity a scales as 1/T 2 . Their ratio is therefore pro- 
portional to temperature, i.e., k/<j ~ T, which is the famous Wiedemann- Franz law. It 
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would be interesting to investigate whether this law also holds here using the holographic 
results for electrical and thermal conductivity. 

As stressed before, the various phases and the corresponding quantum phase tran- 
sition described in section 4 are realized at zero doping. The scale that governs the 
quantum phase transition is an appropriate combination of the parameters g and A ap- 
pearing in our prescription. As a generalization of our holographic description of a Weyl 
semimetal we will extend the prescription to non-zero doping in future work. This is 
interesting for condensed-matter purposes as it corresponds for instance to turning on 
a gate voltage in the material. Holographically this can be done by adding another 
U(l) gauge field to the bulk spacetime under which the fermions are charged, and by 
considering charged Lifshitz black branes [34]. It will be interesting to investigate the 
consequences for the phase diagram described in section 4. 

A topic often discussed in the AdS/CMT literature is the well-known Luttinger 
theorem from Fermi-liquid theory, which states that the total particle density in a system 
is proportional to the total volume in momentum space enclosed by its Fermi surfaces. In 
our case, this can be checked as follows. The system we describe is at zero doping, so the 
total number density of the particles and holes together is zero. The Fermi surfaces for 
particles and holes separately enclose a non-zero volume, but the difference in volume 
enclosed by the Fermi surfaces is zero. Therefore, the Luttinger theorem is trivially 
satisfied in our case. 

Finally, whereas our current prescription takes into account only the leading-order 
contributions in the 1/N expansion, it would be interesting to also consider quantum 
or 0(1/N) corrections. In particular, it is sometimes proposed [35] that the typical 
feature of holographic Fermi liquids, an exponentially suppressed quasi-particle decay 
rate at the Fermi surface, may be a large N remnant, i.e., a consequence of the fact that 
the leading order does not take into account all possible quasi-particle decay processes. 
Therefore, considering also 1/N contributions may resolve this problem and restore 
the conventional u 2 power-law behavior, which is rather robust in Fermi liquids. An 
example of 1/N corrections that are ignored here, are the back-reaction of the massive 
fermions on the gravitational background in the bulk, but one may also think of other 
1/N corrections. 
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A Conventions 



A.l A note on dimensions and units 

In this paper, we work with natural units in which h = c = ks = 1, where ks is 
Boltzmann's constant, and we use dimensionless coordinates r, t, x. As a consequence, 
Newton's constant Gd+i as well as the metric are dimensionless. However, to convert 
the Lifshitz metric (2.1) to standard SI units, we introduce a length scale I and define 
the dimensionful expression 

ds ~ 2 = l 2 ds 2 = - V 2 {f) r —c 2 dP + r -dt 2 , (A.l) 

r 2 V 2 (r) l 2z I 2 

where 

f = lr, x = lx, t = -t, (A.2) 

which has the correct units of meter and second, respectively. Note that I can be 
interpreted as the scale characterizing the size of the Lifshitz spacetime. The temperature 
of the black brane, obtained in (2.2), can be converted to have dimension of Kelvin by 
the rescaling 

T -v B i- < A - 3 > 

and the mass of the fermions in Si-units takes the form 

M 2 = J?- M 2 l d ~ 3 . (A.4) 
Gd+i 

Here, Gd+i is Newton's constant which has dimension m d ■ s -2 • kg -1 . To convert the 
dimensionless frequencies oj and momenta k to SI units, one should rescale 

oj = - uj , k = - . (A. 5) 

Employing the results of this section, the reader can at any time convert to physical 
units. However, for ease of notation, we present our analysis and results in dimensionless 
coordinates. 



A.2 Vielbeins and Dirac matrices 

Here, we collect some formulae related to the spin connection and vielbeins, and state 
our choice of Dirac matrices. The vielbeins are denoted by e- M with fi = 0, . . . , d and 
a = 0, . . . , d, and satisfy 

9p,v = e% rjab e- v , rjab = diag (-1, +1, . . . , +1) . (A.6) 

We will furthermore employ the following notation and relations 

= e%dx^ , n% = dx» , = de^ + n% A , (A.7) 
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where is the torsion-free spin connection. Given the metric (2.1), we can read off 

e-=^- T , e^ = r z Vdt, e l - = rdx\ (A.8) 

TV 

where a = {t, r, i} with i = 2, . . . , d labeling the spatial directions. Using the relation in 
(A. 7) as well as the anti-symmetry of $7, we determine its non- vanishing components as 

ttK = rV d r (r z V) dt , Q,K = rVdx i . (A.9) 

For the spinors, we use the following notation 

* = *t r t 5 ^ = a^ + |(n^, t> = TVv^. (a.io) 

Here, ej^ is the inverse transpose of e- M and the symbol T— is defined in terms of the 
Gamma matrices F- as 

We also perform chiral projections using 

*± = ^ (i±r^-) * , r^± = ±*±. (A.i2) 

A convenient choice of Gamma matrices reads as follows [36]: 
• For odd dimensions d, we can choose 

where the 7-matrices are the <i-dimensional Dirac matrices of the boundary theory. 
We then decompose \£ into chiral components of the bulk spinor 

* - (£) ■ < A -> 

where the components ip± are chiral in the bulk, but they are Dirac spinors on the 
boundary. 

In the case when d is even, a natural choice is given by 

v r - = 1 d+1 , r^ = 7 °, T l - = f-, (A.15) 

where 7 d+1 is the analog of 7 5 in four dimensions and 7 are the gamma matrices 
of the boundary theory. In this case ^± correspond to operators with definite 
chirality on the boundary. 
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• Most of our interest lies in the case d = 4. Using (A. 12), and with a slight abuse 
of notation, we write the four- component Dirac bulk spinor ^ in terms of two- 
component spinors and as 
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(A.16) 



The four-dimensional gamma matrices can be expressed in terms of a- = (1,<?) 
and a- = (— 1, a) with a- being the Pauli matrices in the following way 



pa 







a = {t, i} 



(A.17) 



B Details on the near-boundary asymptotics 
B.l Dirac equation 

The Dirac equation for the fermion field \£ can be obtained from varying the action (2.4) 
with respect to VP, and imposing 5*$>+ = at r = ro (see [14] for more details). We find 
the usual expression (if) — M)^> = 0, which can be written as follows 



rv r r - d r + - r • k + x - r r - Pz (r) - m 



tp(r) = . 



(B.l) 



Here we defined the function p z (r) = r 1 z d r [r z V] + (d — 1)V, and used the notation 
r • k = T- uj + T- k i together with the generalized momenta 

to 



k = (co,ki) , 



uj = — - 



.2-1 



V 



Next, we can simplify equation (B.l) by introducing 

1 



^/ r d-l+zy 



(p(r) , 



which leads to 



rVT r -d r + - T -k - M 
r 



. 



(B.2) 



(B.3) 



(B.4) 



In terms of the chiral components (j)±{r), the Dirac equation (B.4), both in even and 
odd dimensions d, then reads 



<f>±(r) =T- (7- fc) A(TM) <P T (r) 



A(M) = r(rVd r - M) , (B.5) 



where we employed the notation for the T matrices given in appendix A. 2. Finally, we 
can also derive a second-order equation by applying „4(±M)(7 • k) to both sides of the 
Dirac equation in (B.5), which leads to 

k 2 </)±(r) = A(TM)A(±M)<f>±(r) - r 2 V 7 a ( 7 • k)A(±M)<f>±(r) . (B.6) 

k 2 
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Note that for the case of interest in this paper, that is d = 4, half of the components in 
equations (B.5) and (B.6) are trivial, because we can use the Weyl representation of the 
gamma matrices so that <f>± correspond to the upper and lower two components of (p. 



B.2 Asymptotic solutions near the boundary 

We now determine the form of the asymptotic solutions near the boundary at r = oo. 
In contrast to the AdS case with z = 1, for z 7^ 1, the asymptotics are different for 
vanishing and non- vanishing spatial momentum k. We therefore consider both cases in 
turn. 



Asymptotics for k = 

In the case of vanishing spatial momentum, we make a power-law ansatz for the second- 
order differential equation of the chiral components. The asymptotic expansion near the 
boundary then comprises two independent solutions and is given by 

(f> ± = r ±M (1 + r- 2z c± + ...)A± + r^ M ~ z (l + r~ 2z d± + . . .) B± , (B.7) 

where we included subleading corrections in the two separate branches. Here, A± and B± 
are spinors in the Clifford algebra, and the coefficients c± and d± are matrices therein. 
We also observe that (B.7) agrees with the AdS case [36] for z = 1, and we note that a 
subleading term in the A (or B) branch can be more dominant over the leading term in 
the B (or A) branch for certain ranges of M and z. Employing the expansion (B.7) in 
the first order Dirac equation (B.5), we derive a local relation between B T and A± of 
the form 

which fixes half of the integration constants. The coefficients of the subleading terms 
are found as 

2 2 

C± = ; ~ r , d± = — r . (B.9) 

± 2z(±2M-z) ' ± 2z(t2M-3z) v ; 

However, here we need to assume that the subleading terms in the expansion (B.7) are 
still leading over the subleading terms in the Dirac equation (B.5) which appear in the 
function V. A careful analysis translates this assumption into the condition 

z<d-l. (B.10) 

As we can see from (B.8) and (B.9), the above analysis breaks down for M = ±z/2. 
Moreover, in that case the powers in the asymptotic expansion (B.7) can become equal 
and so the spinor coefficients cannot be disentangled. The way to treat the case M = 
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±z/2 is to allow for logarithmic terms in the asymptotic expansion, which lift the de- 
generacy in the two branches A and B. In particular, for M = — z/2 we have 



r 2 A + + r 2 \og(r)B + + . . . , 

r + f (l + c-r~ 2z log(r) + ...^JA- + r~ 3 f (l + ...^jB. 



(B.ll) 



with the coefficients of the form 

B + = i^°A_, B . = i ^A + , C - = Y Z - (R12) 

For M = +z/2, we find 

(f> + = r + ^A + + r~ 3 i B + + ... , 
4>- = r~?A- + r~2 \og(r)B- + . . . , 

together with 



(B.13) 



b ± = r ±M 



B- = -icu 7 °A + , B+ = -^-A-. (B.14) 

2z 

Asymptotics for k 7^ 

For non-zero spatial momenta k and generic values of the Dirac mass M, the asymptotic 
solution to the second-order differential equation (B.6) reads 

(1 + r" 2 C± + ...)A± + r^ M - 1 (l + r~ s± d± + . . .) B± . (B.15) 

Note the important difference compared to (B.7), namely that the exponent of the 
leading order term does not depend on z. The relation between A± and B T can again be 
derived via the first-order Dirac equation (B.5). Together with the subleading coefficient 
c±, their form reads 

The subleading behavior in the B-branch turns out to be more complicated. Both the 
exponent 5± and the coefficient d± depend on the value of z. We have to distinguish 
between three separate cases 



UJ 



/±2M + 1\ _ -. 



1 < z < 3 : S± = z — 1 , d± = — ~~~~r~r 7 " k , 

jfe2 \±2M + zJ 



3<2: * ± = 2 ' rf±= 2(3^2M)- 

As is clear from (B.16), the case of M = ±1/2 is again special and needs to be treated 
separately. Along similar lines as for vanishing momenta, logarithmic terms will appear 
in the expansion. 
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B.3 Hydrodynamic limit 



We can determine the Green's function analytically in the strict hydrodynamic limit 
to — > 0, k — > for arbitrary temperatures T in the region \M\ < z/2. However, one has 
to be careful in defining this limit because the r — > rh and oj — > limits do not commute 
in the differential equation (2.15). Here, we define the limiting procedure as first setting 
k = 0, and then taking oj — > while keeping T finite. Thus, we should use the same 
boundary condition £(r/j) = i as before. The differential equation (2.15) for £(r) then 
becomes 

r 2 Vd r i± + 2Mr£ ± = , (B.18) 

whose solution with the boundary condition ^+(r^) = i reads 

£ ± ( r ) = i e - 2M ^ h ^W . (B.19) 

Substituting then the expression for V(r) given in (2.1) and performing the integral, we 
find 

/r\~ 2M r -, 4M , 

e±W=M-J [1 + F(r)]~3+^ , (B.20) 

where we have left the dimension d and the dynamical exponent z unspecified. The 
self-energy appearing in the Green's function is given by (2.12), and so we obtain 

2M 

G R (0,0) = % -2& Z T- 2 -? , (B.21) 

where we employed (2.2) to related to the temperature T. Note that this expression 
is valid for d + z — 1 ^ 0, which includes the case of interest in this paper, namely d = 4 
and z = 2. 



C Symmetry properties of the Green's function 

In this appendix, we discuss symmetries of the self-energy E(p) defined in (2.12) and of 
the full Green's function (2.11). This analysis utilizes the properties of the differential 
equation (2.15) and boundary condition (2.16), which we recall for convenience 

r 2 Vd r Z± + 2Mrt± = r ^Tk 3 +( r ^ IV ±h)e±, Z±(r h ) = i . (C.l) 

C.l Scale-invariant variables 

We start our discussion with the scaling properties of (C.l). To simplify our notation, 
we define the constant k = d+ ^~ 1 , and we can distinguish three cases which will be 
considered in turn. 
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• We define new variables in the following way 



*=i' t= m- x= k- (c - 2) 



which leads to 



i 



The differential equations (C.l) then takes the form 

a?Vd x Z± + 2Mxi ± = T 1 + (-^ ± l) ^ , e±(^) = i , (C4) 

and we obtain 

£t (x , uJ, T) = fcg 2Af £± (r , w, fc 3 , T) . (C.5) 
Along similar lines, we can define another set of variables as follows 

k=\, T=-, V = 1 t, (C.6) 



leading to 



as well as 



d+z-i — 



y 2 Vd s ( ± +2My( ± = -^-^k+(-^ Y -±k\(l, (±(y h ) = i. (C.8) 

y \y / 

We then find 

U(Vo,k,T) =L>-™Z±(ro,u>,k 3 ,T) . (C.9) 
Finally, we can define variables as 



i 



" = *=(£)"*s, z=^, (CIO) 
where we recall that = (T//e)*. The above choice implies 

v = f^f = ~\ (en) 

and the differential equations become 

z 2 Vd z Z± + 2MzU = ^^Tk+(-^ IV ±k^e±, Z±(z h ) = i. (C.12) 

We then have 

_ 2M 

e±(«o,w,fc) = Q Z £±(r ,a;,fc3,T) . (C.13) 
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C.2 Chirality 



From equation (C.l) we observe that the spin-up and spin-down components £+ and 
are related by a reflection of the momentum 



Z±(r,k 3 ,u;) = ^(r,-k 3 ,Lu). (C.14) 

Furthermore, equations (2.12) and (2.11) which relate £ to the self-energy S(p) and to 
the Green's function £ to Gr do not spoil this symmetry. We thus obtain the following 
relation for the components of the full Green's function 

G%(k,u) = GK-k,u J ). (C.15) 

This symmetry can easily be understood in the relativistic case of z = 1: since ip± 
correspond to left- and right-handed spinors in the boundary theory, the helicity h given 
by h = a ■ k/\a ■ k\ is conserved. In the non-relativistic case z ^ 1, this may be viewed 
as a generalization of chirality. 



C.3 Particle-hole symmetry 

From equation (C.l) we also observe that the components £+ and £_ enjoy the additional 
symmetry 

£±(r,*3,w) = -£(r,-fc 3 ,-w), (C.16) 

where £* denotes the complex conjugate of £. Noting then again that (2.12) and (2.11) 
preserve this property as well, we see that also the full Green's function satisfies 

G%(k,co) = - (G±(-k,-co)Y = - (G*(fc,- W ))* , (C.17) 

where we used the chirality symmetry in the last equality. This symmetry can be seen 
as a particle-hole symmetry since it relates components with frequency oj to components 
with frequency —oj. 



C.4 A symmetry relating M to — M 

In the case when the operator under consideration is right handed, that is C + , the source 
is identified with the boundary value of ip- instead. Hence, the transfer matrix £(r,p) 
defined in (2.7) is replaced with 

C+ = i- — , c --= l i — » c - 18 ) 

where Vi,2 are the up and down components of the spinors defined in (2.5) and (2.6). 
Let us now consider ( + for definiteness, the discussion for is analogous. First, we 
see that £+ = — 1/£+, and we easily verify that the equation for ( + is exactly the same 
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equation (2.15) that £+ satisfies, but with the replacement k 3 — > —k% and M — > —M. In 
addition, the boundary conditions (2.16) for (+ and £+ are the same, and therefore we 
have the following symmetry property: 

(+(r,M,oj,k 3 ) = -^+ 1 {r,M,u,k 3 )=^+{r,-M,u,-k 3 ) • (C.19) 
On the other hand, equation (2.12) is now replaced with 

E(p) = -g lim 2M ((r ,p) , (C.20) 

and so we arrive at 

E(-M,-jfe,w) = -0 2 E _:l (+M,+jfe,u;) . (C.21) 

However, note that strictly speaking this derivation is valid only within the range —z/2< 
M < z/2 for k = and within — ^ < M < ^ for k 0. Outside this region one needs 
to subtract divergent terms to obtain the Green's functions from £ and £, which can be 
shown by observing that the counter terms that are needed in order to renormalize Gr 
also respect this symmetry. 

As a consequence of this symmetry, conditions on the functions s\ t M and s 2 ,m from 
(3.15) can be obtained as 

si-m{u)s 2 ,+m(u) = s 1j+ m{u)s 2 ,-m(u) , ^ 
si,+m(u)si- M (u) - s 2 ,+m(u)s 2 -m(u) = -1 • 

Solving (C.22) in favor of s 2j m(u), we find for non-degenerate values of s^m and s 2 ,m 
that 



s 2 ,m{u) = si,m(h)Ji + i-J 1 i-^r , (C23) 

and by using (C.23) in (3.15) we arrive at 

Gn(k,oj) = -. _ — r . /p r.A\ 

u-\a-k k*' 1 - g k 2M Sl M (u) { 1 + # A /l + r , / n ) 1 j 

Thus, we have reduced the problem to the determination of a single function si : m(u) of 
a single variable u = uj/k z for two values of M. 

C.5 A symmetry of the momentum distribution 

Using the above symmetries of the Green's function in combination with the sum rule, 
we can easily derive a convenient identity for the momentum distribution (4.1). Namely, 
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we have 



AT± + NK = - 

k —k 7r 



J- 



oo 
+00 



ImG^(A;, uj)hf{oj) + ImG^(- k,uj)np{oj) 



(C.25) 



where we used (C.17) in the second line, and the sum rule (2.18) and nir(w)+np(— ui) = 1 
in the third line. 



D Pole structure and sum rules 

In this appendix, we analyze the pole structure of the Green's function, and we show 
analytically for the case of vanishing temperature that the sum rule is satisfied. The 
situation of vanishing momenta can be analyzed for all values of z, however, analytic 
expressions for the Green's function for k 7^ are only available for z = 1. We therefore 
restrict the latter computations to the relativistic case. 



D.l Pole structure 



We now consider the pole structure of the Green's function in the case of vanishing 
momentum and zero temperature. Employing (3.13), we see that there is a trivial pole 
at ijj = which corresponds to the free particle. A non-trivial pole is determined by the 
equation 



= 

Using then polar coordinates oj 
equation 



gc z u 



-i e — (¥+§) 



(D.l) 







with solutions 



sin 



7T 

2 



uj\e td , we find for the imaginary part of (D.l) the 
„ 2M + z it' 



2M 



n + 1 



z 



2M 



irz 



n € Z 



(D.2) 



(D.3) 



If we want to satisfy the Kramers-Kronig relations, that is the requirement of no poles 
in the upper half-plane, we have to make sure that there is no solution to (D.3) in the 
range (0, +7r) for all values of n, that is 



9n i (0, +tt) 



(D.4) 



By noting that < z — 2M < 2z for \M\ < |, we see that the only possible solution 



in the forbidden range is given by n = — 1 with B-\ 



All the other solutions are 
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already out of the forbidden range, but yield poles on the non-principal sheets of the 
complex plane. For example, for n = — 2 the pole resides in 6-2 € (vr, +00) depending 
on the value of M, while for n = it is in the range #0 £ (0, —00). 

Let us now turn to the real part of (D.l). Employing again polar coordinates and 
using (D.3), we arrive at the equation 



= 1 - gc z 



, 2M 1 , 

M" (-1) 



which is solved by 



\gcz 



sign(£(C 2 ) = (-1) 



(D.5) 



(D.6) 



It is now clear that in order to avoid the pole Q-\ = | on the principal sheet, we have 
to require 



gc z > . 



(D.7) 



Furthermore, all other possible poles of the Green's function (specified by (D.3) and 
(D.6)) with even n are either on the lower-half of the principal sheet or on non-principal 
sheets of the complex plane. Therefore, for the choice of sign (D.7) the Kramers-Kronig 
relations are satisfied. 

For the case of vanishing temperature and z = 1, analytic results also for k / 
are available. As it has been explained in [14], the analytic continuation of the Green's 
function for non-vanishing momenta is given by replacing ui — > \/uj 2 — k 2 . This implies 
that there are poles at ui = ±|fe|, with the pole in the lower half-plane corresponding 
to the retarded Green's function. Furthermore, for the choice of sign (D.7), there are 
no additional poles in the upper half-plane of the principal sheet. Therefore, we have 
shown that in the relativistic case the Kramers-Kronig relations are satisfied also for 
non-vanishing momenta. 



D.2 Sum rule for k = 

The Green's function in the case of vanishing temperature and momenta can be found 
in equation (3.13). The spectral density (2.17) then reads 



p(jfe = 0,w) = ^-lmTi\G R (k = 0,u)] = — Im 

2"7T 7T 



Introducing the following constants 

B = gc z cos 
we express the spectral density as 







Kt 









lo — gc z uj z e 



gc z sin 



2M ■ 1 M , 1 \ 







Kt 









p(k = 0,w) = + 



C LO~ 



(l>- Bu,-) + (Coj—\ 



2 • 



(D.8) 



(D.9) 



(D.10) 
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Next, to compute the sum rule we recall the symmetry p(k, —uj) = p(k, +oj). We then 
have 

"+oo r+oo 

du p{k = 0, oj) = 2 / da; p(k '■ = 0, uj) 
Jo 

2 z 

~ n 2M - z 

2M 



f 

j — i 



arctan 



B '2mB 2 + C 2 ^ =0 ° 



(D.ll) 



J u>=0 



sign (g c z ) . 



2M - z 2M - z 

Therefore, for the choice of sign (3.3) we indeed find that the sum rule (2.18) for vanishing 
momenta and vanishing temperature is satisfied 

»+oo 

dcup(k = Q,u) = 1 . (D.12) 



D.3 Sum rule for k ^ 

In the case of z = 1 and vanishing temperature, the expression for the Green's function 
was given in equation (3.1) together with (3.4). To compute the sum rule in this situation, 
we write 



/+oo 
-oo 



dui p(k, ui) 



r+oo 

Jo ^ 



du Im - 



-2oj 



27f "p 2 (l - g Cl e-^^+lyM-l^j ' 

with p as in (3.4). Next, we take (3.4) into account to obtain 



(D.13) 



I r+oo 

H / dwlm 

* U\k\ 



-2oj 



(oj 2 - k 2 ) 2 (l - gey e -^(^+D( W 2 _ fc 2) 



2M-1 



(D.14) 



where we note that the first term on the right-hand side vanishes due to taking the 
imaginary part. Performing then the change of variables x = co 2 — k 2 , we find 



/+oo _^ roo 
du)p(k,u) = / dx Im 
-oo ^ Jo 



x 2-g Cl x 2M+l e -i*(M+l) 



After a further change of variables y = x 2 , we arrive at 



/_" 



+oo _^ 2 r°° 

du p(k, oj) = / dy Im 

n Jo 



y-g Cl y 2M e -iAM+\) 



(D.15) 



(D.16) 



which is a special case of the computation in appendix D.2. We therefore find 

/+oo 
dup(k,uj) = 1 . (D.17) 
-oo 
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E WKB calculation 



In this appendix, we present details of the calculation leading to equation (4.22). We 
consider the case M > 0, T = and without loss of generality make the assumptions 
to > 0, k = (0,0,^) and > 0. Furthermore, we consider the spin-up component of 
the Green's function, the other cases can be worked out using the symmetries outlined 
in appendix C. 

The parameter which we take large in the WKB approximation is k = k/y/uj, and 
the fluctuation equations of u± are given by (B.6) which yield 



= d?u± + 



H =r- ± k) 

r 2 k 2 \ r J 



2 

r r 2 k 2 
= d 2 r u± + fi(r)d r u± + ff{r)u± , 



d r u± + 



fc 2 



r 3 k 2 \ r J 



M(M± 1) 



u± 
(E.l) 



where k had been defined in (B.2). This expression can be put in the Schrodinger form 
by the transformation 



u± = h±(r)y± , h±(r) = exp (~\ J ^/fW^ • 



(E.2) 



A suitable parametrization of the extra spatial coordinate in the limit k 3> 1 is given in 
terms of x = r/^fuj. In this variable, the transformation (E.2) reads as follows 



h±(x) -- 

Then, the Schrodinger equations become 

1 d 2 



Vxk±l 



c 3/2 



k 2 dx 2 

with the effective Schrodinger potentials 



y±(x) - V±(x)y±(x) = , 



V±{x) 



1 



x^ 



k 2 x 6 



+ 



a±M)(l±M) : 3 1 



k 2 x 2 



+ - 



A(xk±l) 2 xk(xk±l) 



(E.3) 



(E.4) 



(E.5) 



In addition to solving the above Schrodinger equation for y + analytically using the WKB 
approximation as described below, we have also solved it numerically. The numerical 
results are presented in figure 18, which we discuss as we proceed in our computation. 

We now focus on the plus component of the fluctuations. In the region close to the 
horizon at x -C 1, and in the limit k S> 1, the relevant piece of the potential for y + 
reduces to 



V. 



1 



+ 



x 



4 p x 6 ' 



(E.6) 
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(a) Double-logarithmic plot of \V+ 
point is approximately x tp = 1/10. 



The turning 



(b) Double- logarithmic plot of Rey + . The imagi- 
nary part is similar but with a minus sign. Left of 
the turning point there is oscillatory behavior, the 
other regions exhibit exponential behavior. 




0.10 



(c) Near-horizon oscillatory behavior. Until the turning 
point, there is good agreement between the numerical (blue) 
and WKB (red) result. After that, the WKB solution is no 
longer valid. 




0.11 



(d) Matching procedure near the classical turning point. The 
dashed red curves are the WKB results, the solid blue curve 
is the numerical solution, and the solid green curve is the 
Airy function solution needed in the matching formulae which 
agrees with the numerics around the turning point. 
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Rey + 

10 39 r 




0.1 0.5 1.0 5.0 10.0 50.0 100.0 



(e) Double-logarithmic plot of Rey+ in the classically disallowed re- 
gion. The solid blue curve is the numerical solution, the dashed lines 
are the analytic solution (E.12) for exponents a = n/2 (red) and 
a = 5/3 (yellow). The former agrees with the numerical result, the 
latter deviates from it. For large x, the WKB results are not valid 
anymore but are taken over by the Bessel-function solutions plotted on 
the far right, which overlap with the WKB solution for smaller x. 



Rey+ 




10 1 00 1 000 10 4 10 5 1 6 



(f) Double-logarithmic plot of the near-boundary region, where the 
Bessel-function solution is valid. The blue curve is the numerical result, 
the dashed curves are the analytic result (E.16) for exponents a — tv/2 
(red) a = 5/3 (yellow). 

Figure 18: Illustration of the numerical WKB calculation for k — 10 and M = 1/4. 
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from which we see that there is a turning point located at approximately x tp = 1/k. In 
the region x < xt p the WKB wavefunction is given by 

v .(x) = — e +ik ^ tp d W- v +W + ^ e ~ fk ^ tp d V-^+W (E 7) 

with V+ obtained from (E.6). The solution for x < xt p is compared with the numerics 
in figure 18(c). This solution should be connected to the solution near the horizon at 
x — > 0, which yields the condition that A2 = 0. We will not need to fix the coefficient 
A\ because, as shown below, it will cancel out in the final answer. The next step is to 
continue the solution from the region x < xt p to the region x > x tp by using the well- 
known WKB connection formula, see for example [25, 26]. This procedure is illustrated 
in figure 18(d). One important subtlety here is that in the final answer we will need 
second-order corrections to the WKB formula as well. Noting that the parameter that 
plays the role of h here is 1/k, we find, 



A ie - 1 * 

y+{x) 



V+(x) 



O-(X) 



(E.8) 

where the second-order WKB coefficients can be found as (see for example [37]) 

^ = ±rJ\YM*W\. (E . 9) 

This is our WKB solution in the classically disallowed region x tp <C x < 00, where we 
can now approximate the potential by 6 

v + ^\ + m±^ + ^L-^. ( E.io) 

Uj rv Jb Ju i\i Ju i\i 

To the desired order in k, (E.9) then yields 

±a± = J_ + M{M + 1)*> MJ_1_ {En 
± 4k 2 x 4k 2 8k 3 4k ' 

and with this approximation we find from (E.8) the following solution, which is illustrated 
in figure 18(e), 



y + (x) = A 1 xe ! 4 ye * e k k + {x) + -e + * e k K-(x) J , (E.12) 
where we employed 



6 Note that we keep the leading 1/k 2 corrections in order to be consistent with keeping the a± 
corrections in (E.8). 
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and have denned the functions 



K± 



M(M + l)x M - 
1± 2k ± W l ° g{kx) 



M(M + l)x 2 
1 Ak 2 + 



Ax 2 k 2 



^ 4& 2 4A; 3 x 



(E.14) 



Note that here we only show terms up to the desired 1 jk 2 order. The terms in the first 
brackets come from the corrections to the exponential term in (E.8), the second terms 
come from the expansion of the V + ^ 4 and the last ones from the a± corrections in (E.8). 

The next step is to connect this solution to the region near the boundary at x — > oo 
and read off the Green's function from the coefficients of the near-boundary expansion. 
Note that the WKB approximation in the near-boundary limit x — > oo fails, because in 
this region we have an inverse-square potential which does not satisfy the WKB condition 
\Vl(x)/V2(x)\ < 1. Thus, one has to solve the Schrodinger equation with the potential 
(E.10) exactly and connect it to the solution (E.12) in the overlapping region k/x ^> 1. 
In the latter, we can approximate the potential as 



V+(x) 



1 + M(M 



1) 



x^ 



x 2 k 2 



(E.15) 



and the exact solution to the Schrodinger equation with this potential is found in terms 
of the Bessel functions as 



y+(x) = ^fx 



+ B 2 Ji. 



■M 



(E.16) 



where are coefficients so far undetermined. This solution is also shown in figure 
18(e). In order to connect it to (E.12), we use the well-known asymptotic formula for 
the Bessel functions for large values of the variable k/x and we find 



y+(k/x) 



x 



V2nk 
from which we obtain 

B^ 



B 1 (el + e"^ M e-i) + B 2 (e§ - e^ M e "l) 



+ 



(E.17) 



Aie~ i ^V2nk 



AttM —k' 



2cos(vrM) V 2 

We shall not need the expression for B 2 , and the constants e± are defined as follows 



(E.18) 



M T 1 

e± = 1 ± ^73- log k =F 



2k 2 



Ak 2 ' 



(E.19) 



On the other hand, using equation (E.3) together with the series expansion of the Bessel 
functions near x — > 00 given in (E.16), we find that 



u+(x) ^A u + x M + ... , 



(E.20) 
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where the superscript u denotes that we consider the spin-up component of the spinor 
A+. This is in accord with the asymptotic solution given in appendix B.2 with 



Al = B 



l-M 2 M+\ 

n\-M) ' 



(E.21) 



where B\ is determined by (E.18). 

The next step in the calculation of the Green's function is to obtain the WKB solution 
for the u- component in (E.l). One way would be to apply the same steps as above in 
the WKB calculation, however this process is more difficult than the u + case, since the 
potential for u_ shown in (E.5) has extra poles at x = 1/k. Fortunately, u_ is related 
to u+ by the first order differential equation (B.5), and so we can use the result for u+ 
to obtain u_ directly. Equation (B.5), specified to the spin-up component and in the 
rescaled variables x and k, reads 



U- 



k + 



T (x 2 d x - Mx) u + , 



which for the solution (E.8) becomes 

2 f h '+( X ) 



U. 



k + 



1 V'ix) 



h+{x) 4 V(x) 



Mx 



u + + A(x) 



(E.22) 



(E.23) 



where the function A{x) is given by 



A(x) 



IX 



2 h+ A ie - in / 4 



k+l V l J\x) 



[kV_?(x) 
(kvj(x) 



1 + 
1 + 



k ) k 

cr_(x)\ cr'_(x) 



k 



k 



(E.24) 



On the other hand, very close to the boundary, we can approximate the Schrodinger 
potential as 



V-(x) 



1 + M(M 



I) 



2 k 2 



(E.25) 



and again, using this potential in the Schrodinger- like equation for y_ , the exact solution 
reads 



y-(x) = jx 



(E.26) 



with C\ : 2 some coefficients to be determined. Similarly as above, expanding the Bessel 
functions in this solution for large x/k, and using U- = h-(x)y-(x) with (E.3), we find 



u-(x) ^A a _x- M + ... , 



(E.27) 
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with 



A u = d = . (E.28) 

Next, we employ the asymptotic formula for the Bessel functions in (E.26) for large k/x 
to match to the WKB result (E.23). In this matching, we also need to use (E.24), (E.12) 
and the relation u±(x) = h±(x)y±(x) with (E.3). The result is 

2cos(7rM) \2 / 
and we shall not need the expression for C2. The constants j± read 

7 ± = l±^log* T j|. (E.30) 

Finally, as can be inferred for instance from (2.7), the spin-up eigenvalue of the matrix 
£ is proportional to the ratio of u_ and u+. Thus, using its definition from (2.12), the 
spin-up component of the self-energy is proportional to the ratio of the constants ^4" 
and A\ . To be more precise, 

£+ = -igu M S- = -igT 2M k 2M T - M \ ^ . (E.31) 

If we then use (E.18), (E.21), (E.28) and (E.29), we obtain the final result for the spin-up 
component of the self-energy as it was given in (4.21). 

Before closing this appendix, let us remark the following. As explained at the begin- 
ning of section 4.3, the numerical prefactors both of the terms in e± and 7± and the one 
in the exponent, i.e., the number a in the end result ImS~ exp(— ak 2 /uj), are difficult 
to determine analytically. This is because the integral in(E.8) leading to n± cannot be 
done analytically. In the calculation presented here, we have approximated the poten- 
tial V+ by taking into account only the three most relevant terms in the region x > xt p , 
which leads to a = 2 as in (E.12). This is reasonable away from the turning point, but 
close to the lower integration limit x tp the term has a considerable contribution not 
taken into account. Furthermore, the turning point itself is determined by setting to 
zero an approximation to the effective potential. To do better, we should use the full 
potential everywhere, and expand the integrand as a series in 1/k. However, this series is 
not uniformly convergent because the integration limit xt p depends on 1/k. In practice 
we can terminate the series at some point, which leads to a different approximation of 
the numerical prefactors. For example, if the full effective potential is taken into account 
and the integrand y/V+ is expanded in 1/k, the next contribution of order k 2 is —k 2 /6, 
which results in a = 5/3, and thus in a decay rate of F ~ exp(— 5fc 2 /3a;). Furthermore, 
we can do the same for the integral in the a± from (E.9). Also, in the prefactor V + ^ 4 
the full potential can be taken into account. Then, the 0(1/ k 2 ) terms in j± and e± also 
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Figure 19: Plot of i dty/V+(t) as a function of k, which appears in the exponents of (E.8) 

that determine the quasi-particle decay rate. For large k, the k 2 term in the exponent is dominant 
and therefore the curves converge to the numerical prefactor of this term. In this way, the correct 
factor can be determined numerically. Here, we have taken M = +1/4 and x — 20 as the upper 
integration limit, which is an intermediate value of x. The value of the integral is independent 
of M for k — > oo. 

The solid green curve is the result when the two most relevant terms in the potential are taken 
into account, leading to a — 2 (dashed green line). The solid blue curve is the resulting expression 
when the integration over the full potential is done analytically before terminating the series in 
k, leading to the a = 5/3 prefactor of the k 2 term (the dashed blue line is the number 5/6). The 
red dots show the value of the integral when it is performed numerically, and the dashed red line 
is the number 7r/4. The former converge to the latter curve, which is evidence that the value of 
the integral is approximately ir/Ak 2 + . . .. Thus, the exponent that determines the quasi-particle 
decay rate in the final WKB result will be exp(— irk 2 /2uj). 



receive extra corrections. Still, we are left with an approximation of the exact numbers, 
and it is not clear whether there occur cancellations between terms later on in the 1/k 
expansion. 

The integral over y/V+ can however be performed numerically. Figure 19 compares 
the numerical and analytic values of the number a/2 that appears in the exponent. The 
numerical result is approximately 7r/4, whereas our improved approximated analytic 
calculation yields 5/6. When we use a = n/2 in the exponent, the resulting self-energy 
agrees more accurately with the numerical results, as is shown in figures 15 and 18. 
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